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ABSTRACT 

The main uncertainty in current determinations of the power spectrum normalization, erg, from abundances of X-ray luminous 
galaxy clusters arises from the calibration of the mass-temperature relation. We use our weak lensing mass determinations of 
30 clusters from the hitherto largest sample of clusters with lensing masses, combined with X-ray temperature data from the 
literature, to calibrate the normalization of this relation at a temperature of 8 keV, Msooc.gkev - (8.7 + 1.6) h^^ 10^* Mq. This 
normalization is consistent with previous lensing-based results based on smaller cluster samples, and with some predictions from 
numerical simulations, but higher than most normalizations based on X-ray derived cluster masses. Assuming the theoretically 
expected slope a = 3/2 of the mass-temperature relation, we derive erg - 0.88 ± 0.09 for a spatially-flat ACDM universe with 
Q„, - 0.3. The main systematic errors on the lensing masses result from extrapolating the cluster masses beyond the field-of-view 
used for the gravitational lensing measurements, and from the separation of cluster/background galaxies, contributing each with 
a scatter of 20%. Taking this into account, there is still significant intrinsic scatter in the mass-temperature relation indicating 
that this relation may not be very tight, at least at the high mass end. Furthermore, we find that dynamically relaxed clusters are 
75 + 40% hotter than non-relaxed clusters. 

Subject headings: Cosmology: observations — cosmological parameters — dark matter — gravitational lens- 
ing — Galaxies: clusters 



1. INTRODUCTION 

The abundance of massive clusters of galaxies provides sen- 
sitive constraints on the cosmological parameters that gov- 
ern structure growth in the universe. However, a prerequisite 
for this is reliable mass measurements for large samples of 
clusters with well-understood selection criteria. Cluster mass 
measurements used have traditionally come from virial anal- 
ysis of the velocity dispersion measurements of cluster galax- 
ies (e.g., Frenk et al. 1990; Carlberg et al. 1997; Borgani et al. 
1999), or X-ray temperature measurements of the hot intra- 
cluster gas under the assumption t hat the gas is in hydrostatic 
equilibrium (for a review see iRosati. Borgani & NormanI 
l200a) . Satellite observatories such as ROSAT, ASCA, XMM- 
Newton and Chandra have made increasingly accurate X- 
ray temperature measurements of clusters, and have produced 
well-defined cluster samples of sufficient size to accurately 
measure the X-ray temperature and luminosity functions (e.g. 
lHenryl2()04l : lBohringer et alJ2002h . However, the relation be- 
tween cluster mass and X-ray temperature and luminosity, re- 
spectively, must be determined to convert these into a reliable 
cluster mass function. 

X-ray luminosities are available for large samples of clus- 
ters, but the X-ray luminosity is highly sensitive to the com- 
plex physics of cluster cores, making it challenging to re- 
late to cluster mass. Measuring X-ray temperatures is ob- 
servationally much more demanding, but the X-ray tem- 
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perature is mainly determined by gravitational processes 
and is hence more directly related to cluster mass than X- 
ray luminosity. Both from simulations and observations 
the intrinsic scatter in mass around the mass-temperature 
relation is thus found to be much smal ler (AM/M =a 
0.15, [E vrard. Metzler, & NavaiTO 1996; Bor gani et al]l2004l: 
[Sanderson et alJI2003t IVikhUnin et al.l i2006) than the scatter 
in mass around the mass-lu minosity relation (AM/M ^ 0.4, 
iReiprich & Bohringeiil2002h . 

Two different routes have been followed for determining the 
mass-temperature relation. Most studies have used a small 
sample (up to about a dozen) of supposedly well understood 
clusters for which the assumptions underlying the mass de- 
termination should be fulfilled to a high degree. The main 
concern for this approach is that the selected clusters may not 
be representative of the whole cluster population, and there- 
fore the derived mass-temperature relation may only apply 
to a subset of clusters. Alternatively, the mass-temperature 
relation may be determined from a large sample of more ob- 
jectively selected clusters. This is more fruitful when com- 
paring such a locally determined mass-temperature relation 
to a sample of high-redshift clusters where the data quality 
does not allow a similar selection of the "most suitable" clus- 
ters. Also, mass-temperature relations derived from simula- 
tions are usually based on a large range of simulated clusters 
with no pre-selection. Hence it is most appropriate to com- 
pare observationally obtained mass-temperature relations de- 
termined from all available clusters to the relations from sim- 
ulations. On the other hand, for some of the clusters in such 
a sample the hydrostatic assumption may be invalid, making 
X-ray based mass determinations unreliable for a subset of 
the clusters. A larger scatter (which may not be symmetric) 
around the mean mass-temperature relation may be expected, 
when such clusters are included. 
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There are still poorly understood systematic uncertain- 
ties associated with establishing the mass-temperature re- 
lation. The normalization of the mass-temperature rela- 
tion based on cluster masses deter mined fro m X-ray data 
jFinoguenov, Reiprich, & Bohringer 2001; iVikhlinin et alJ 
I2OO6; Arnaud. Pointecout eau, & Pratt 2005). tend to differ 
significantly between studies, and from the expectations based 
on numerical simulations (e.g., Evrard et al. 1996; Eke et al. 
1998; Pen 1998; Borgani et al. 2004). The determination of 
this normalization is currently the dominating source of dis- 
crepancies between the reported values for the power spec- 
trum normalization on the scale of galaxy clusters, erg, derived 
from the observed cluster temperature function (Huterer & 
White 2002; Seljak 2002; Pierpaoli et al. 2003; Henry 2004). 
Observations using X-ray based mass determinations have 
traditionally favored low normalizations (and hence low val- 
ues of cTg), while simulations have favored somewhat higher 
normalizations. 

Gravitational lensing provides an opportunity to measure 
cluster masses without invoking the assumption of hydrostatic 
equilibrium in the hot intra-cluster gas implicit in the X-ray 
based mass determinations. Also, in this case the measure- 
ment of cluster mass is truly independent of the X-ray tem- 
perature measurement. Hjorth, Oukbir & van Kampen (1998) 
used weak gravitational mass measurements for eight clusters 
drawn from the literature to find a relation between mass over 
cluster-centric radius, and temperature. They determined a 
normalization of this relation consistent with the value pre- 
dicted by Evrard et al. (1996), but with a preference for some- 
what higher cluster masses (if the redshift scaling of equa- 
tion|5]is assumed, see below). However, Smith et al. (2005; 
hereafter S05) determined a mass-temperature relation with a 
normalization significantly lower than indicated by the Hjorth 
et al. (1998) study. S05 based their results on a sample of 
10 clusters with weak lensing masses and temperatures deter- 
mined from Chandra data. 

Here, we present a new weak gravitational lensing-based 
measurement of the normalization of the mass-temperature 
relation. The main improvements with respect to the work 
of Hjorth et al. (1998) and S05 is that we use a significantly 
larger cluster sample which represents a significant fraction 
of all the clusters in an even larger sample with well-defined 
objective selection criteria (Dahle et al. 2002; H. Dahle 2006, 
in preparation). An additional improvement over the work 
of Hjorth et al. (1998) is that the weak lensing analysis has 
been performed in a consistent way for all clusters, using 
the same shear estimator and making the same assumptions 
about e.g., the typical redshift of the lensed galaxy population 
and the degree of contamination by cluster galaxies. We note 
that the early data set of clusters with published weak lensing 
masses used by Hjorth et al. (1998) is biased at some level to- 
wards systems that were observed because of "extreme" prop- 
erties, such as being the hottest or most X-ray luminous sys- 
tem known at the time, or having a large number of strongly 
gravitationally lensed arcs. Furthermore, we note that our 
gravitational lensing measurements are made at larger radii 
than probed by S05, requiring smaller extrapolations to esti- 
mate the mass within e.g., the virial radii of the clusters. 

The data set used for the analysis is described in § |2] our 
results for the mass-temperature relation and erg are presented 
in § [3] and our results are compared to other work and the 
implications discussed in § H] 

Except when specifically noted otherwise (for easy compar- 
ison to previous results using different cosmologies), we as- 



sume a spatially-flat cosmology with a cosmological constant 
(Q,„ = 0.3, Oa = 0.7), and the Hubble parameter is given by 
Ho = lOO/ikms-'Mpc^^ 



2. DATA SET 
2.1. Weak lensing data 

Our weak lensing data set is a sample of 30 clusters (see Ta- 
ble [TJ, of which 28 were included in the weak lensing cluster 
sample of Dahle et al. (2002). Data for two additional clusters 
come from a recent extension of this data set (H. Dahle 2006, 
in preparation). The clusters targeted for these weak lensing 
studies were generally selected to lie above an X-ray luminos- 
ity limit Lz,o.i-2.4keV > 6 X 10'*"' ergs s"' (this luminosity limit 
is for our chosen cosmology with h = 0.7) and within a red- 
shift range 0.15 < < 0.35. The observed clusters were 
selected from the X-ray luminous cluster samples of Briel 
& Henry (1993) and Ebeling et al. (1996;1998;2000). The 
cluster samples of the first two of these papers are based on 
correlating an optically selected cluster sample (Abell 1958; 
Abell, Corwin, & Olowin 1989) with X-ray sources from the 
ROSAT All-Sky Survey (RASS; Trumper et al. 1993), while 
the two last papers contain X-ray flux limited cluster catalogs, 
also based on RASS. Of the total sample of 30 clusters, three 
(A959, A 1722, and A 1995) are drawn from the Briel & Henry 
(1993) sample and two (A209 and A2104) are drawn from the 
XBAC sample of Ebeling et al. (1996). Of the remaining 25 
clusters, 22 are included in the X-ray brightest cluster sample 
(BCS) of Ebeling et al. (1998), while three (A611, A1576, 
and Zw3146) come from its low-flux extension (eBCS; Ebel- 
ing et al. 2000). Of the BCS and eBCS clusters in our sam- 
ple, 24 objects are included in a volume-limited sample of 
35 clusters selected from the BCS and eBCS samples (Dahle 
2006). Hence, while our sample is not strictly physically well- 
defined (in the sense that the availability of an X-ray temper- 
ature measurement is one of the defining selection criteria), it 
still has significant overlap with a well-defined cluster sam- 
ple. In a recent paper, Stanek et al. (2006) discuss how a 
significant scatter around the mean mass-luminosity relation 
may cause a significant Malmquist bias in X-ray flux-limited 
cluster samples, causing high-mass, low flux clusters to drop 
out at high redshifts. This would result in a bias in the mass- 
luminosity (or mass-temperature) relation derived based on 
such a sample. We note, however, that although the RASS- 
based samples from which our cluster sample is drawn are 
flux-limited, the cluster sample discussed here quite closely 
approximates a volume-limited sample, and we therefore ex- 
pect any such bias to be negligible. 

The observations were made with the 8192^ UH8K mosaic 
CCD camera and the 2048^ Tek CCD camera at the 2.24m 
University of Hawaii Telescope and with the 2048^ ALFOSC 
CCD camera at the 2.56m Nordic Optical Telescope. All clus- 
ters were imaged in both the /- and V-band, with typical total 
exposure times of 3.5h in each passband for the UH8K data 
and ~ 1 .5h for the data obtained with the more sensitive 2048^ 
detectors. The seeing was in the range 0'.'6 < FWHM < I'.'l 
for all the imaging data used for the weak lensing analy- 
sis. The median seeing was 0'.'82 in the /-band and 0'.'9 in 
the V-band. This gave typically ~ 25 usable background 
galaxies per square arcminute, or a "figure of merit" value 
of 2 Q^/dQ. ^ 1.5 X 10-**deg-2, as defined by Kaiser (2000). 
As noted below, the background galaxies were selected based 
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on signal to noise ratio rather than magnitude, with limits cor- 
responding to 21 < mi < 24.5 and 22 < my < 25.5 for point 
sources. The observations and data reduction of the data set 
used for the weak lensing mass measurements are described 
in detail by Dahle et al. (2002). 

Major efforts are being made to improve the methods for 
the estimation of weak gravitational lensing, particularly in 
connection with ongoing and future studies of "cosmic shear" 
based on wide-field optical surveys. The requirements for the 
precision of shear estimates in these surveys are substantially 
more stringent than for weak lensing observations of massive 
clusters, given the significantly weaker lensing effects mea- 
sured in random fields. 

In this work, we have used the shear estimator of Kaiser 
(2000), which was "blind-tested" (along with several other 
shear estimators) by Heymans et al. (2005), using simulated 
lensing data. The shear estimator of Kaiser (2000) is more 
mathematically rigorous than the currently most widely used 
shear estimator (Kaiser, Squires, & Broadhurst 1995), but 
it displays a significant non-linear response to shear, unlike 
most other shear estimators. If we correct our shear values 
using a second order polynomial based on the test results of 
Heymans et al. (2005), we find that most cluster masses stay 
within +/ - 15% of the mass calculated based on uncorrected 
shear values. Furthermore, the change in average cluster mass 
is < 2%, i.e., there is very little systematic shift in mass. In the 
end, we chose not to apply this correction, since it would, in 
a few cases, require extrapolations outside the range of shear 
values over which the shear estimator has been tested. For 
more details about the practical implementation of this shear 
estimator, see Dahle et al. (2002). 

To convert the measurements of weak gravitational shear 
into actual cluster masses, the distances to the background 
galaxies need to be known. The background galaxy redshifts 
were estimated from spectroscopic and photometric redshifts 
in the Hubble Deep Field (for details, see Dahle et al. 2002). 
For our data set and chosen cosmological model, the average 
value of the ratio between the lens-source and observer-source 
angular diameter distances, /3 = Di JDs, is well approximated 
by the relation {/?) = 1.37z^,-2.00zd + l-01 within the redshift 
range of our cluster sample. This then provides an effective 
critical surface density for lensing {l,ait = (c^ /4-nG)(Di(J3}y^ ; 
where Di is the angular diameter distance to the cluster), 
which is used for deriving cluster masses from the shear es- 
timates. The quoted value of (J3} corresponds to the value at 
large cluster radii; at smaller radii a correction term has to be 
employed to account for contamination by cluster galaxies, as 
discussed below and illustrated in Figure [T] 

The observable galaxy shape distortions caused by gravi- 
tational lensing provide a measurement of the reduced tan- 
gential shear, gj - jt !()■ - k), where yj is the tangential 
component of the shear and k is the convergence. We fit an 
NFW-type mass density profile. 



ScPciz) 



(1) 



(r/r,)(l + r/r,)2 

(Navarro, Frenk, & White 1997), to the observed reduced 
shear profile gjir) of each cluster. Here, pdz) is the critical 
density of the universe at the redshift of the cluster, and 
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(2) 



3 ln(l + C2oo) - C2oo/(l + C2oo) 
We assumed a concentration parameter C200 - Cvii-/1.194 - 



4.9/(1 -I- z), corresponding to the median halo concentration 
predicted by Bullock et al. (2001) for a Mvk ^ 8 x IO'^Mq 
cluster from simulations of dark matter halos in a ACDM uni- 
verse. Here, C200 = rioadrs, and Cvii- = r^ulrs, where r2ooc is 
defined as the radius within which the average mass density is 
200 times the critical density pdz), and ryir is the virial radius 
of the cluster 

The lensing properties of the NFW model have been cal- 
culated by Bartelmann (1996) and Wright & Brainerd (2000). 
From our fit, we calculated Msoor, the mass enclosed by the ra- 
dius rsooc. The mass estimates are listed in Table[T] The shear 
measurements used for the fit were made at clustercentric 
radii 50" < r < 180" for the clusters that were observed with 
2048^ CCD cameras and 150" < r < 550" for the clusters that 
were observed with the UH8K camera. By comparison, we 
find rgooc values typically in the range 300" < rgooc < 600" 
for the clusters we study here. In many cases, we need to 
extrapolate the NFW profile out to rgoor (in Table [T] we list 
the ratio of the outermost radii of our shear measurements, 
rfit, to rgooc, and note that rgooc - 0.66r2ooe for our chosen 
NFW model). In this extrapolation, we assume the median 
NFW concentration parameter given above. Hence, any in- 
trinsic scatter in Cyir will introduce an extra uncertainty in the 
cluster mass estimates. If we assume a random scatter around 
the mean value of Cvir at the level (a Icr A(logCvir) ~ 0.18) 
predicted by Bullock et al. (2001), we find a corresponding 
scatter in the mass estimates of 20% for our data set. This 
additional scatter is not included in the uncertainties of the 
listed mass measurements in Table [T] but is considered fur- 
ther in Section 3.2. 

The measured gravitational lensing signal is sensitive to the 
two-dimensional surface mass distribution, including mass 
associated with the cluster outside rgooc, and random struc- 
tures seen in projection along the line of sight (Metzler, White 
& Loken 2001; Hoekstra 2001; Clowe, De Lucia & King 
2004; de Putter & White 2005). This will inti'oduce addi- 
tional uncertainty (and potentially a net bias) to any lensing- 
based estimates of the cluster mass contained within a 3D vol- 
ume. Studies based on simulated clusters (e.g., Clowe et al. 
2004) indicate that the net bias is no more than a few percent 
when 3D cluster masses are estimated by fitting observations 
of grir) to predictions from theoretical models of the mass 
distribution, such as the NFW model. However, the scatter in 
the mass estimates from projection effects amount to a weak 
lensing mass dispersion of ~ 15-25% for massive galaxy 
clusters, which should be added to the observational uncer- 
tainties of the lensing mass estimates. In this paper, we have 
assumed a lensing mass dispersion of 0.26 resulting from pro- 
jection effects, corresponding to the value estimated by Met- 
zler, White & Loken (2001) from their N-body simulations. 
Although these authors considered a somewhat different mass 
estimator, more recent estimates indicate a similar mass dis- 
persion for the NFW profile fitting method that we have used. 
This additional mass uncertainty has been added in quadrature 
to the uncertainties of Mgooc values listed in Table [T] 

The absence of reliable information about the individual 
redshifts of the faint galaxies used for the weak lensing mea- 
surement will inevitably result in some degree of confusion 
between lensed background galaxies and unlensed cluster 
galaxies. The magnitude of this effect will depend on the 
projected number density of cluster galaxies, and should thus 
have a strong dependency on cluster radius. Hence, a radi- 
ally dependent correction factor was applied to the shear mea- 
surements to correct for contamination from cluster galax- 
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Fig. 1. — Cluster galaxy contamination in the faint galaxy catalogs as a 
function of distance from the cluster center The solid line represents an 
average of 6 clusters at an average redshift {z) = 0.31, while the dashed line 
represents an average of 5 clusters at {z) = 0.23. 

ies in the faint galaxy catalogs that were used to measure 
the gravitational shear (these catalogs included all galaxies in 
the cluster fields that were detected at a signal-to-noise ratio 
6 < S /N < 100, with no additional selection based on e.g., 
galaxy color). The magnitude of this correction was estimated 
from the radial dependence of the average faint galaxy density 
in two "stacks" of clusters observed with the UH8K camera, 
one at z ~ 0.30 and the other at z ~ 0.23, assuming that the 
contamination is negligible at the edge of the UH8K fields, 
> 1.5/i ' Mpc from the cluster center. The estimated degree 
of contamination is shown in Figure [T] 

Given the difference in cluster redshift, the similarity of the 
two curves in Figure [T] may be somewhat surprising, as one 
would naively expect the more distant clusters to display a 
significantly lower surface density of cluster galaxies. How- 
ever, there are several competing effects that affect the ob- 
served galaxy density at a fixed angular radius: Firstly, if all 
cluster galaxies were detectable regardless of cluster redshift, 
the change in apparent image scale with redshift should in- 
crease the surface density by a factor given by the square of 
the ratio of the angular diameter distances. On the other hand, 
a fixed angular radius would correspond to a larger physi- 
cal cluster-centric radius (and hence lower galaxy density in 
physical units) at the larger redshift, the difference depend- 
ing on the slope of the radial galaxy density profile. At the 
radii probed in this study, both the radial surface mass density 
profile and the number density profile of bright cluster galax- 
ies follow approximately the power law behavior of a singu- 
lar isothermal sphere (cr(r) oc r", with a - -1). Hence, the 
physical number density (in galaxies/Mpc^) at a fixed angular 
radius should decrease as the inverse of the ratio of angular 
diameter distances. In addition, the faintest galaxies drop be- 
low the detection limit at higher redshift, the eff'ect depend- 
ing on the slope of the cluster luminosity function around 
Mr X -15. Assuming a Schechter (1976) luminosity func- 
tion with a faint-end slope a* - -1.25 (typical of rich clus- 
ters) and M*{z - 0.23) = -21.65, the luminosity function 
can be integrated down to the detection limit (corresponding 
to Mr ^ -15.0 and Mr ^ -15.7 at z = 0.23 and z = 0.30, 
respectively), to estimate the fraction of cluster galaxies that 
drop out at the higher redshift (~ 30%). Finally, a redshift- 
dependence given by M^(z) = M|(0) + 5 log(l + z) was as- 
sumed to account for galaxy evolution in the clusters. The 
combination of all these effects would predict a surface den- 



sity of cluster galaxies which is 7% less sA z - 0.3, com- 
pared to z = 0.23, for a fixed cluster richness. Even this 
small difference would be erased by a slight decrease in the 
assumed values of the slopes a* and a. A faint-end slope of 
the luminosity function of a* = -1.1 would be sufficient to 
remove the predicted difference in galaxy surface density at 
the two different redshifts. Based on Hubble Space Telescope 
(HST) WFPC2 imaging of the galaxy cluster A2218 (which 
is similar to the clusters studied here in terms of optical rich- 
ness, lensing mass and X-ray properties), Pracy et al. (2004) 
find that the cluster core shows a relative depletion of dwarf 
galaxies, leading to a radial profile of faint galaxies which is 
significantly shallower than the SIS prediction. For the "inter- 
mediate" dwarf population (-18 < Mpeoew < -15, similar to 
the range in absolute magnitude of cluster galaxies in our faint 
galaxy catalogs), these authors find a radial distribution with 
a slope a - -0.63 + 0.09. Assuming a similar slope for the 
radial distribution of the faint cluster galaxies in our catalogs 
would also remove the predicted difference between the two 
curves in Figure[T] Based on this figure, and the above discus- 
sion, we conclude that we are probably justified in ignoring 
the redshift-dependence in our cluster galaxy contamination 
correction. 

The level of cluster galaxy contamination for individual 
clusters will generally differ from the mean level calculated 
above, as there will be significant cluster-to-cluster variations 
in the abundances of cluster dwarf galaxies. Based on the 
sample of clusters observed with the UH8K camera, the scat- 
ter in dwarf galaxy richness was estimated to be ~ 50% (this 
estimate also includes variations in the field galaxy density 
caused by uncorrelated large-scale structures along the line of 
sight, and hence the true scatter in dwarf galaxy richness of 
the clusters is somewhat overestimated). By employing the 
mean contamination correction calculated above rather than 
an estimate appropriate for each cluster, we introduce an ad- 
ditional scatter of up to 20% in our mass estimates. This ad- 
ditional scatter is not included in the uncertainties of the tabu- 
lated mass measurements in Table[Tl but is considered further 
in Section 3.2. 

Eight of the clusters in our sample were also included in 
the combined strong and weak gravitational lensing study of 
S05, based on observations of a sample of 10 X-ray lumi- 
nous galaxy clusters at z ~ 0.2 using HST WFPC2. These 
authors estimated the projected cluster mass within a cluster- 
centric radius of 250/;"' kpc, assuming an Einstein-de Sitter 
(Q,„ = 1, Qa = 0) cosmology. Figure |2] shows a compar- 
ison of the mass values listed by S05 with our cluster mass 
estimates, using the best-fit NEW model to derive projected 
cluster masses, assuming the same cosmology as S05. These 
authors assumed a spatially constant contamination of 20% 
cluster galaxies in their background galaxy catalogs at radii 
< 2', while we find an average contamination of 30% for our 
data at these radii. Hence, for the plot in Figure|2]we have ad- 
justed our radially dependent contamination correction such 
that the average contamination at small radii is consistent with 
that assumed by S05. We find that our cluster mass estimates 
are generally consistent with those of S05, although with a 
tendency for higher masses (by about 30%). 

2.2. X-ray data 

For clusters in the weak lensing data set, we compiled 
a list of corresponding X-ray temperatures from the litera- 
ture. For many of these clusters, their global temperature, 
or even a temperature map, has been determined using data 
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Fig. 2. — Comparison with estimates of the projected mass within 250/j"' 
kpc published by SOS. The lensing masses in this plot are given for the cos- 
mological model assumed by SOS, and our estimate of the degree of galaxy 
cluster contamination was also adjusted to match the contamination assumed 
by SOS (see text for details). 

from Chandra and/or XMM-Newton. However, these tem- 
peratures constitute a rather heterogeneous sample, for which 
the systematics are not well established. Consequently, X- 
ray temperatures were primarily dra wn from the samples of 
lota & Mitsudal (|2004 . lAUenl (lioOOl) and [White (2000), each 
providing a homogeneous measure of the global cluster tem- 
perature (i.e., temperature measured within a cluster-centric 
distance close to rgooc) for a large fraction of the clusters in 
the weak lensing sample. All these authors derived tempera- 
tures based on analysis of ASCA spectra. For two clusters not 
in either of the samples mentioned above we extracted pub- 
lished temperatures from other sour ces. 

Specifically, from the works of Ota & Mitsuda' (2004") and 
[Mien (2000) we extracted temperatures estimated by these 
authors by fitting an isothermal plasma model with the Galac- 
tic absorbing column densit y as a free parameter. The tem- 
peratures taken from Whit3 (l2000h were derived by fitting an 
isothermal plasma model with the nominal Galactic absorb- 
ing column density fixed, but since only energies above 1 keV 
are used in the I White! (|2000) spectral fitting the fixed col- 
umn density should not intro duce systematic efects relative 
to the tempe ratures from the lOta & Mitsudal (l2004t) and the 
lAUenl (12000 ) samples. For the remaining two clusters we ex- 
tracted published temperatures obtained in a similar way (see 
TablelD. 

For the clusters in two or more samples, their derived tem- 
peratures agree within the uncertainties, and for any clus- 
ter the derived temperatures differ by less than 20% be- 
tween samples. Also, the mean of temperat ure dif f erence s 
between any two of the samples by Ota & Mitsuda (2004), 
lAUenl (i2000) , and IWhitei (.2000) is less than 3%, indicating 
the low level of systematic temperature variance between dif- 
ferent analyses. 

Although the isothermal plasma model has proven too sim- 
plistic for nearby clusters, the global cluster temperature is 
straightforward to derive from observations as well as simu- 
lations, enabling a rather direct comparison between obser- 



vations and theory. Furthermore, for the majority of distant 
{z > 0.5) clusters only global isothermal temperatures can be 
obtained in the foreseeable future. Hence, we refrain from go- 
ing into the detailed spatial and spectral modeling of the intra- 
cluster gas. The effects of cluster dynamics, "cooling cores", 
non-sphericity etc. generally aff ects the global temperatures 
only at the 10%-20% level (e.a:..rEvrard. Metzler. & Navarrol 
fT996HWhitai2000l:ISmith et al.ll2005l) . 

3. RESULTS 
From the virial relation 

M oc {v^)r oc Tr, (3) 

between cluster mass, M, inside radius r, galaxy velocity, v, 
and gas temperature, T, combined with the definition of mass 
within an over-density of 500 times the critical density 

MsQOc °^ rl^cPc(z) °^ 4)0c'E'^(^)Pc(0) t>c rl^^E^iz), (4) 

where the term E^(z) - pc{z)lpc(Q) describes the evolution of 
the over-density for a given cosmology, the mass-temperature 
relation is obtained as 

E{z)Mioo, oc T^/\ (5) 

In this study, we have measured Mgooe, relating to Mvir 
through Mvii - 1 .65M500f for the average cluster redshift with 
our chosen NFW model. We take into account that the slope 
may deviate from the simple theoretical expectation a - 3/2 
and normalize the relation at 8 keV since our sample is dom- 
inated by massive clusters. Hence, mass-temperature rela- 
tions were obtained by fit ting the data in Table [T] us ing the 
BCES(X2|Xi) estimator of lAlcritas & Bershadvl (1 19961) to the 
following parameterization of the mass-temperature relation 

£(z)M5ooc = M5ooc,8kev(r/8keV)«. (6) 

The redshift-dependent factor, E{z), contained in eq|6]must 
be calculated individually for each cluster, as this would oth- 
erwise produce an artificial 15% variation in mass over the 
redshift range spanned by our cluster sample. In effect, 
the normalization of the relation refers to the present epoch 
(z = 0). 

The fitting procedure of lAkritas & Bershadvl (Il996h takes 
uncertainties in temperatures as well as in weak lensing 
masses into account, and makes no assumptions about the 
intrinsic scatter of both quantities. Results from fitting sub- 
samples as well as the full sample are presented in Table |2] 
and Figure [3] 

For the full data set (for those clusters with temperature 
from more than one sample the tem p eratur e was taken i n 
prioritized order from lOta & Mitsudal (|2004 . lAllenI (l2000l) . 
an d lWhita(l2000l) 'l we find the following normalization of the 
mass-temperature relation at 8 keV Msooc.skev = (8.7 + 1.6) x 
IO'^^/i^'Mq and a slope of cc = 0.49+0.80. It is evident that the 
slope of the mass-temperature relation is not well-determined 
since our data only span a modest range at the high mass/high 
temperature end of the cluster distribution. In fact, it is not 
obvious that there is a tight mass-temperature relation at the 
high temperature end probed here. 

3.1. Normalization of the mass-temperature relation and cr^ 

The concentration of clusters around 6-8 keV enables 
a robust measurement of the normalization of the mass- 
temperature relation at the high-mass end. Even though the 
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TABLE 1 

Weak lensing masses and X-ray temperatures 



Cluster 



2d" 



rut /'•500c 



kT (Ota & Mitsuda) 
(keV) 



kT (Allen) 
(keV) 



kT (White) 
(keV) 



kT (other)" 
(keV) 



7.54 
8.79 
8.67_ 
25.27 
3.83 



A68 0.255 21.37+ 

A115 0.197 2.42 ' 

A209 0.206 

A267 0.230 

A520 0.203 

A586 0.171 

A611 0.288 

A665 0.182 

A697 0.282 

A773 0.217 

A959 0.285 

A963 0.206 

A1576 0.299 

A1682 0.226 

A1722 0.325 

A1758N 0.280 

A1763 0.228 

A1835 0.253 

A1914 0.171 

A1995 0.320 

A2104 0.153 

A2111 0.229 

A2204 0.152 7.86_ 

A2219 0.228 4.72^^ 

A2261 0.224 11.52f 

MS1455+22 0.258 3.21^^ 

RXJ1532.9+3021 0.345 13.86^ 

RXJ1720.1+2638 0.164 3.28 ' 

RXJ2129.6+0005 0.235 8.35 

Zw3146 0.291 7.57 



5.40+g 
12.18!; 
13 09+*"' 
9.33+3-5" 
4.42 
8.62 

2.70!?i^ 
20.37+|3^ 
4 90+2 ''^ 

8 42+'' '" 

2 62+2 03 

23 69+*-*^ 
Ml 



14.14 
3.84 



0.30 
0.50 
1.62 
1.32 
1.20 
0.21 
0.59 
0.36 
0.39 
0.30 
1.45 
1.51 
1.58 
0.55 
2.58 
0.33 
0.45 
0.39 
0.44 
1.23 
0.23 
0.47 
0.27 
0.45 
0.32 
0.58 
0.44 
0.39 
0.38 
0.47 



6.93+"*3 



<: r 1 +0.44 
■'■■^ -0.41 



6 96+" '''' 



6 96+"-28 

8 19+0-62 
°-'^-0.60 
Q 07+0.70 
=■"'-0.66 
=; 24+0-89 
-0.73 

6.83!OJ[ 



6.42+0-63 

6.88!0.«6 
8.11+0?:* 



7.42+ 



9 06+'" 

' -1.32 

7.66+0-''5 
6 94+8:"^^ 
6.68+8| 



6.56+ 



4.91+ 



5.72+ 



7 94+0-96 

7 02+^-^^ 
' -"^-0.80 



o n+0.62 



8.29+ 



6 B+O ''^ 

°-^-'-0.30 



7 33+0-35 

-^-0.30 



6 23+0-30 

°-^-'-0.28 
9 46+0-63 
^■^°-0.57 

°-°^-0.46 
:i+0.27 
0.25 



4.33^ 



6 06+0'''' 

-Q52 
6.85!; 

7 70+0.41 

8.63+!'-«« 



6.08! 
6.57+ 



7.24+ 



7 30+0 *6 

-0.38 
7 00+0.49 

10 53+0 5' 

J jy+I.OT 

Q 1 9+8:I§ 
'■'^-0.46 



6.99+ 



Q+0.24 
'-0.23 

9 52+0-" 

7 49+0.57 



4.83 



-0-43 
+0.22 
-0.21 



7.10+ 



5.60+ 



6.80!0-3S 



" — See Ebeling et al. (1996; 1998) for referenc es to redshift measur ements.' 
nication); kT value for RXJ1720. 1+2638 from lMazzotta et alj <200lft . 



- kT value for A209 from P.B. Marty (private commu- 



1 0.0 F 




4 6 8 10 12 14 

kT, (keV) 



Fig. 3. — Weak lensing mass estimates and X-ray temperatures with 
BCES(X2|Xi) regression lines for different subsamples. The full line is the 
fit to the full sample, defined as described in the text; the dotted line is the 
fit to t he iOta & Mitsudai 12004) sample; the dashed lin e is the fit to t he lAHer] 
120001) sample; and the dot-dashed fine is the fit to the I White! <2000l) sample. 

slope of the mass-temperature relatio n varies substantially 
between the three X-ray sub-samples (Ota & Mitsuda' l2"004t 
[Alien .2000; White 2000 ) the best fit normalizations agree 
within their statistical uncertainty. We note that for all fits. 



the four different regressions of lAkritas & Bershadvl (Il996h 
all result in normalizations within 20%. 

The strongest constraints on the mass- 
temperature relation normalization are obtained 
by taking advantage of previous studies of mas- 
sive clusters ^noguenov, R eiprich. & Bohringe^ 

200lt lAUen. Schmidt & FabiaS lloOlt 

Arnaud. Pointecouteau. & Pratt! |20()5Ik showing that the 
mass-temperature relation slope is close to a = 3/2 
as expected from simple gravitational collapse models 
dKaiserl 11986 ). Hence, in order to express the normal- 
ization of the mass-temperature relation in terms of the 

characteristic temperature Tj, - 8 k eV(1.6 5M5Q0c.8kev) " 
dPierpaoli. Borgani. Scott. & White! l2003h we assume 
(X = 3/2 (with a representative un c ertain ty of 10%, 
e.g. iFinoguenov. Reiprich. & Bohringed 1200 Ih . As is 
custom for quoting values we adopt the redshift 

dependence factor F(z) = (AcE^)' ^'^ [l - 2^]"^^^ 
(iPierpaoU. Borgani. Scott. & Whitd l2003h where is the 
mean overdensity inside the virial radius in units of the 
critical density at the relevant redshift {F(z) and E(z) differs 
by 7% at z = 0.23). 

We find Ti, - 1.28 + 0.20 for our full sample and re- 
sults from calculations of based on various subsamples 
are listed in Table |3] From the erg - relation plotted by 
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iPierpaoli. Borgani. Scott. & Whitd (l2003h in their Figure 2, 
we find erg = 0.88 + 0.09, based on our full sample. We 
note that this relation is valid only for an intrinsic scatter in 
temperature of < 10% around the mean mass-temperature re- 
lation. A larger intrinsic scatter will imply a lower value of 
erg. We provide our constraints on the intrinsic scatter below. 

3.2. Scatter in the mass-temperature relation 
The squared scatter in lensing mass, M500C, around the best 



fit, {o-'"'f 

surement error, the squared intrinsic scatter, and the squared 
systematic errors (cT"^f = {a'lff + {cr'^f + (crj)^. The 
main systematic errors in the lensing mass (see Section 2.1) 
arise from extrapolating the assumed NFW mass profile out 
to rsQQc (due to cluster-to-cluster variations in the assumed 
concentration parameter c,,„ ) and from the separation of clus- 
ter/background galaxies (due to cluster-to-cluster richness 
variations). Each of these introduces a scatter of 20% in the 
lensing mass, hence o"^' = 0.28. For the full sample we find 
cr'^ - 0.94 which is larger than expected from the mean lens- 
ing mass error, cr^'' - 0.52 and the systematic errors, indicat- 
ing either a sizable intrinsic scatter in mass or that the mea- 
surement/systematic errors are severely under-estimated. 

Accounting for errors in both mass and temperature, we 
find an intrinsic scatter in of 0.25 ^^ jg. There is a 70% 
probability that the scatter in temperature is larger than 10%, 
favoring somewhat lower values of erg than quoted above. 
However, most of the scatter is caused by the low mass clus- 
ters. 

3.3. Relaxed vs. non-relaxed clusters 

We looked into whether relaxed clusters and non-relaxed 
clusters have the same normalization of the mass-temperature 
relation. Our "relaxed" cluster sample consists of A586, 
A963, A1835, A1995, A2204, A2261, RXJ1720, and 
RXJ1532. These are clusters with "spherical" optical and X- 
ray morphology, and no known cluster-scale dynamic distur- 
bances. For the relaxed clusters we find a normalization of the 
mass-temperature relation of Msooc.gkev.reiax = (17.3 + 3.7) x 
lO''*/! 'M0 while the normalization for the non-relaxed clus- 
ters is M5()0c,8kev,nonreiax = (7.6 + 1 .5) X 10''*/i"'Mo (scc Fig- 
ure |4||. The higher normalization of relaxed clusters is sup- 
ported by the fact that the mean mass of relaxed clusters is a 
factor 1 .5 larger than the mean mass of "non-relaxed" clus- 
ters, although the relaxed and the non-relaxed clusters span 
roughly the same temperature range 

The scatter in mass for the relaxed sample (cr^ - 
similar to the scatter for the non-relaxed sample (cr'^' 
The mean error for both samples is cr^'' = 0.57. Either relaxed 
clusters spread as much around their mass-temperature rela- 
tion as clusters in general, or we have used a poor definition 
of "relaxed" clusters. However, the fact that the normalization 
of the mass-temperature relation for relaxed clusters is higher 
than for non-relaxed clusters indicates that there is a physi- 
cal difference between the two sub-samples. From the present 
study, it thus seems that relaxed clusters do not form a tighter 
mass-temperature relation than clusters in general. 

For a given mass, non-relaxed clusters are found to be 
~ 75 + 40% hotter than relaxed clusters. Since we con- 
sider global, isothermal temperatures, the presence of "cool- 
ing cores" in relaxed clusters will result in a lower global tem- 
perature than the virial tem perature. However , this eff'ect is 
at the 10%-20% level (e.g.. Ismith et"ani2003aL lO^Hara et all 



<(- 



— ) ), is the sum of the squared mea- 



- 0.77) is 
= 0.88). 




Fig. 4. — Weak lensing mass estimates and X-ray temperatures with 
BCES(X2|Xi) regressions for "relaxed" clusters (filled symbols, full line) and 
"non-relaxed" clusters (open symbols, dashed line). 



TABLE 2 

Best-fit mass-temperiCTUre relation (ARBrraARY slope) 



-'WjoOcSkeV 


a 


Sample 


(lo'^/r'Mo) 






12.0 ±2.4 


1.30 ±0.97 Ota<< 


i Mitsuda (2004) 


9.5 ± 1.9 


0.87 ± 0.78 


Allen (2000) 


8.1 ± 1.5 


0.11 ±0.98 


White (2000) 


8.7 ±1.6 


0.49 ± 0.80 


All 


17.3 ± 3.7 


2.69 ± 1.30 


"Relaxed" 


7.6 ± 1.5 


0.29 ± 0.76 


"Non-relaxed" 



TABLE 3 

Best-fit mass-temperature relation (fixed 
slope) 



7"* 



Sample 



1.04 ±0.18 3/2 Ota & Mitsuda (2004) 

1.22 ±0.21 3/2 Allen (2000) 

1.35 ±0.21 3/2 White (2000) 

1.28 ±0.20 3/2 All 

0.82 ±0.14 3/2 "Relaxed" 

1.42 ±0.22 3/2 "Non-relaxed" 

Note. — The uncertainty in from a has 
been taken to be 10%. 



120061) so this cannot alone explain the temperature differ- 
ence between relaxed and non-relaxed clusters. Based on the 
mass-temperature relation from 10 clusters (3 of which are 
considered relaxed), S05 also find that non-relaxed clusters 
are hotter than relaxed clusters. An objective classification of 
the degree of relaxation for a sizeable cluster sample is re- 
quired for further quantifying the size of this effect. 

4. DISCUSSION 

Based on the hitherto largest sample of X-ray luminous 
clusters with measured lensing masses, we derive a normal- 
ization of the mass-temperature relation at the high mass end, 
A^500c,8kev = (8.7 + 1.6) /i"' IO^^Mq. This value is higher 
than the lensing based mass-temperature normalization of 
S05, based on a smaller cluster sample, but is consistent with 
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TABLE 4 

Normalizations of the mass-temperature relation 



Method 


(z) 


-'^^SOOc.SkeV 


Slope 


Reference 


X-rays 


0.09 


6.00 ± 0.35 


Fitted 


Amaud. Pointecouteau, & Pratt (2005) 




0.09 


6.07 ± 0.46 


Fitted 


VikhUnin et al. (2006) 


Lensing 


0.23 
0.23 


6.65 ± 0.52 
8.67 ± 1.57 


Fixed 
Fitted 


S05 
This study 


Simulations 


0.04 
0.00 


8.09 ± 0.48 
6.82 ± 0.76 


Fixed 
Fitted 


Evrard, Metzler. & Navarro (1996") 
Borgani et al. (2004) 



Note. — Included here are some recent studies with 
a quoted normalization of the M500C - T-relation, scaled 
to MsoocSkcVi the mass contained within rsoor of a 
cluster with kT = 8keV, with eiTors propagated, and as- 
suming h = I. The listed normalization has been scaled 
to a common redshift of 2 = 0. The normalization of 
S05, which was determined within a fixed physical radius 
of 250/i" ' kpc in an Einstein-de Sitter universe, has been 
scaled to M500C f™ our adopted cosmology, making the 
same assumptions about an NFW-type mass profile as we 
have made for our own data. For each study, we indicate 
whether the slope of the relation was calculated from a fit 
to the data, or whether the slope was fixed at the theoreti- 
cally expected value, a= 3/2. 



this within Icr errors; see Table |4] Mass-temperature rela- 
tions with masses determined from X-ray data tend to have 
a lower normalization than lensing based relations, and they 
are only marginally consistent with our nor malization. This 
is also the case for the two recent studies of IVikhlinin et alj 
(120061) and lArnaud. Pointecouteau. & Pratt! (l2005h based on 
smaller sample s of lower mass (and hence cooler) clus- 
ters. [Vikhlinin etd] (I2OO6I) measured cluster masses inside 
rsooe from X-ray observations of a sample of 13 low red- 
shift clusters with a median te r nperat ure of 5.0 keV while 
lArnaud. Pointecouteau. & Pratt! (!2005h determined the nor- 
malization from X-ray derived masses of 10 nearby clusters 
with a mean temperature of 4.8 keV. The two studies agree 
on the same normalization, higher than previous X-ray mass 
based studies, but there still seems to be a ~ 20% discrep- 
ancy between X-ray and lensing derived mass-temperature 
relations. 

We note, however, that the lensing based and X-ray based 
normalizations are made at different redshifts, and that this 
discrepancy would vanish if the redshift-dependence pre- 
dicted by the self-similar collapse model in equation |5] were 
neglected. Given the heterogeneous nature of these data 
sets, any claim of significant departures from self-similarity 
would be premature, but this clearly provides an interesting 
avenue for future research, involving even larger cluster sam- 
ples spanning a wider interval in redshift. 

We confirm the result of Smith et al. (2005) that non- 
relaxed clusters are on average significantly hotter than re- 
laxed clusters. This is qualitatively consistent with N- 
body/hydrodynamical cluster simulations which show that 
major mergers can temporarily boost the X-ray luminosities 
and temperatures w ell above their equilibrium values (e.g. 
!Randall et al]|2002i) . 

In contrast to several previous (mainly X-ray mass based) 
published mass-temperature relations, the normalization de- 
rived in this study is in good agreement with the normalization 
derived from numerical simulations. However, the accuracy 
of the normalization is not good enough to discriminate be- 
tween simulations including different physical processes. Our 



results show that X-ray based measurements of the cluster 
abundances, after reducing the major systematic uncertainties 
associated with the mass-temperature normalization, give an 
amplitude of mass fluctuations on cluster scales that is con- 
sistent with other methods. This lends additional support to 
the "concordance model" cosmology, and lends credence to 
the basic assumptions of Gaussian density fluctuations. Our 
determination of erg - 0.88 ± 0.09 is higher than most erg 
determi nations fr om cluster data (for a compilation of these, 
see e.g.. !Henryl2004.) . However, our finding is consistent with 
the value derived from weak gravitational lensing in the com- 
bined Deep and Wide CFHT Legacy Survey (erg = 0.86+0.05; 
Semboloni et al. 2005) based on the halo model of density 
fluctuations (Smith et al. 2003b). It is also consistent with 
the CMB+2dFGRS+Ly„ forest result (erg = 0.84 ± 0.04) 
of Spergel et al. (2003), with the joint CMB + weak lensing 
analysis of Contaldi, Hoekstra, & Lewis ( 2003), which gave 
erg = 0.89 + 0.05, and with CMB analyses ( !Bond et al.!!2"005l) 
yielding erg a 0.9. However, the more recent 3-year WMAP 
results (Spergel et al. 2006) give a significantly lower value 
of erg, and also a preference for a value of Q„, lower than 
0.3. Also, results from the recent 100 square degree weak 
lensing survey (Benjamin et al. 2007) favor a lower value of 
erg = 0.74 + 0.06 for Q,„ - 0.3. We note that our quoted 
value of erg is based on the assumption that the intrinsic scat- 
ter about the mass-temperature relation is < 10%, and that 
our erg estimate will be biased high if the true scatter signifi- 
cantly exceeds this value (Pier paoli. Borgani. Scott. & Whitd 

llool . 

The limiting factor of our measurement of the normaliza- 
tion of the mass-temperature relation is the magnitude of the 
measurement errors (dominating the systematic errors, esti- 
mated to be ~ 30%). In order for the mass-temperature re- 
lation to be a competitive route for constraining cosmologi- 
cal parameters and to discriminate between simulations with 
different input physics, the normalization must be measured 
to better than ~ 10% accuracy. However, there are good 
prospects for improving on these results in the near future. 
Firstly, the superior spectro-imaging capabilities of Chandra 
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and XMM-Newton will allow the construction of large, ho- 
mogeneous cluster temperature samples. A comparison to 
tailored simulations with realistic physics, analyzed in the 
same way as observations, will advance our understanding of 
systematics and the link between the mass-temperature re- 
lation and structure formation (C.B. Hededal et al. 2007, in 
preparation). Secondly, more accurate weak lensing-based 
mass measurements of a larger sample of clusters are feasible 
as large mosaic CCD cameras that can probe intermediate- 
redshift clusters beyond their virial radii are now common, 
and the cluster sample could easily be doubled from a similar 
survey in the Southern celestial hemisphere. 

Finally, we note that a more direct measurement of erg from 
weak lensing by clusters is possible, provided that weak lens- 
ing mass estimates are available for a large, well-defined, 
volume-limited cluster sample. This makes it feasible to 
calculate the cluster mass function directly from the lensing 
masses, rather than indirectly via the X-ray temperature func- 
tion (Dahle 2006). Since mass estimates based on baryonic 



tracers of the total cluster mass only enters indirectly as a se- 
lection criterion (e.g., clusters selected based on X-ray lumi- 
nosity above a certain threshold), the method is less suscep- 
tible to systematic and random errors, as it does not require 
an accurate characterization of the scatter around the mean 
mass-temperature relation. 
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